A first-principles study of the electronic structure and stability of Be(BH4)2 
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Alanates and boranates are studied intensively because of their potential use as hydrogen storage 
materials. In this paper we present a first-principles study of the electronic structure and the 
energetics of beryllium boranate, Be(BH4)2. From total energy calculations we show that - in 
contrast to the other boranates and alanates - hydrogen desorption directly to the elements is 
likely, and is at least competitive with desorption to the elemental hydride (BeH2). The formation 
enthalpy of Be(BH[4)2 is only —0.12 eV/H2 (at T — OK). This low value can be rationalized by 
the participation of all atoms in the covalent bonding, in contrast to the ionic bonding observed in 
other boranates. From calculations of thermodynamic properties at finite temperature we estimate 
a decomposition temperature of 162 K at a pressure of 1 bar. 
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I. INTRODUCTION 

In the last decade the environmental importance of re- 
ducing the CO2 exhaust has been widely accepted. The 
use of hydrogen based fuel cells is an important contri- 
bution to achieve this reduction. One major obstacle for 
this use is the development of a method for hydrogen 
storage with a high gravimetric and volumetric hydrogen 
densityji 

One way of storing hydrogen is in a (complex) metal 
hydride. The ideal hydrogen storage material should 
have the highest possible gravimetric hydrogen density. 
This obviously requires the use of lightweight materials. 
Moreover, the formation energy of the hydride has to be 
such that it is stable at atmospheric conditions, yet it has 
to decompose at a moderate temperature to release the 
hydrogen. A further important point is that the reac- 
tions involved in hydrogen de/ absorption must have fast 
kinetics. 

Over the last decade alanates and boranates have 
been studied extensively as potential hydrogen storage 
materialsjii^ Alanates and boranates consist of a lattice 
of metal cations and (Alll4)^ or (BIl4)^ complex an- 
ions, respectively. Generally these materials decompose 
by heating via intermediate complex hydrides into bulk 
metals, elemental hydrides and hydrogen gas. In the 
last few years the attention has gradually shifted from 
alanates towards boranates, because of the high gravi- 
metric hydrogen density in the latter. Many boranates 
turn out to be too stable, however. 

In principle a large variety of boranates can be syn- 
thesized by changing the metal cations, which can be 
used to tune the formation energy;^, So far most effort 
has been devoted to the alkali boranates,i^i^i^i^i2iiiS and 
more recently to mixtures of alkali boranates^^^ and 
to the alkaline earth boranateS i^'^'^'^d^d^d'^ In order to 



understand the chemical trends we have recently de- 
veloped a simple model for the formation energies of 
these compounds. This model demonstrates that these 
boranates are ionic compounds (in the sense discussed 
above) and that the difference in their formation ener- 
gies can be understood on the basis of the electrostatic 
(Madelung) lattice energy. The basic stability of the 
(BIl4)~ cation is not affected by substituting one alkali 
or alkaline earth cation by another. 

The stability of (BH4)~ may be changed by adding an 
element that competes with boron in binding with hy- 
drogen. To investigate this possibility we study beryl- 
lium boranate, Jie(B}i^)2r^ in this paper. Establish- 
ing the electronic structure and themodynamic stability 
of Be(BIl4)2 will assist in understanding the chemical 
and physical trends in alkali, alkaline earth alanates and 
boranates i2£ 

We present a first principles study of the elec- 
tronic structure and the thermodynamic properties of 
Bc(BII^)2. The electronic structure in relation to the 
crystal structure is used to analyze the bonding in 
Be(BII^)2. We calculate total energies and phonon fre- 
quencies of all compounds involved in possible formation 
reactions of Be(BH4)2. From these data we obtain the 
thermodynamic properties at finite temperature. 



II. COMPUTATIONAL METHODS 

First-principles calculations are carried out within the 
density functional theory (DFT) approach, applying a 
generalized gradient approximation (GGA) for the ex- 
change correlation functional. We use a plane wave 
basis set and the projector augmented wave (PAW) 
methodf^Sii^ as implemented in the Vienna Ab initio 
Simulation Package (VASP),— '^^1^^ and apply non-linear 
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core corrections <21 

Brillouin zone integrations are performed with a tetra- 
hedron method^^ for calculating total energies. A Gaus- 
sian smearing method is used for calculating densities of 
states, with a smearing parameter of 0.1 eV. The k-point 
meshes are such that total energies are converged within 
0.1 meV per formula unit. The total energies used in 
the calculations of the reaction enthalpies are calculated 
with a high plane wave kinetic energy cutoff of 700 eV. 
By varying the computational parameters, in particular 
by trying different PAW potentials,— we estimate that 
reaction enthalpies are converged on a scale of 5 meV. 

The atomic positions and lattice parameters are re- 
laxed using a conjugate gradient algorithm for a range of 
fixed volumes. The total energy versus volume curve ob- 
tained this way is fitted with a Murnaghan's equation of 
state expression, which yields the ground state volume, 
the bulk modulus, and its pressure derivative.— At the 
ground state volume we relaxed the atomic positions and 
lattice parameters to obtain the ground state structure. 
This procedure is followed for all compounds mentioned 
in this paper. 

To calculate the zero point energies (ZPE) and phonon 
densities of state we need the phonon frequencies of all 
these compounds. Vibrational frequencies are obtained 
from the dynamical matrix, whose matrix elements (the 
force constants) are calculated using a finite difference 
method^i The force constants are calculated from dis- 
placements of 0.005 A in two opposite directions for each 
atomic degree of freedom. For both bulk beryllium and 
beryllium hydride 2x2x2 supercells give converged 
ZPEs. One does not need a supercell to calculate the 
phonon frequencies of Be(BH4)2, since the unit cell of 
Be(BH4)2 is sufficiently large. For boron we use the fre- 
quencies that have been reported earlier 




FIG. 1; Schematic bonding scheme and labeling of the atoms 
within a polymer chain in 3e{B}i^)2- The three dimensional 
structure is given in Fig. [2l The numbers indicate optimized 
bond lengths in A. 




III. CRYSTAL STRUCTURE 



FIG. 2: (Color online) Crystal structure of Be(BH4 



Be(BH4)2 can be synthesized by the reaction of lithium 
boranate and beryllium chloride. Its crystal structure 
consists of helical polymers of alternating beryllium and 
boron atoms (Bf,) that are connected via pairs of hydro- 
gen atoms (Hf,)»^ The polymer building block is shown 
schematically in Fig. [1] A further boron atom (B^) is 
attached to each beryllium atom, again via a pair of hy- 
drogen atoms (He) and this B^ atom also binds two "dan- 
gling" hydrogen atoms (Hd). The polymers are packed 
in the crystal structure as shown in Figure [H On the 
basis of this structure one may expect a strong bonding 
between the atoms in one polymer chain, and a much 
weaker bonding between the polymer chains. The lat- 
ter is reflected in the low melting point of Be(BH_j)2 of 
125°C. 

We relaxed the crystal structure of Be(BH4)2 as de- 
scribed in the previous section, including the cell vol- 
ume, lattice parameters and atomic positions, while keep- 
ing the experimental space-group and Wyckoff positions. 



Except for the cell volume, the parameters compare well 
to the experimental values.— The calculated cell volume 
is 15% larger than the experimental cell. This indicates 
that the binding between the polymer chains is indeed 
weak and of van der Waals type. It is well-known that, 
using the common functionals, DFT fails to capture van 
der Waals bonding. However, the total energy difference 
between the experimental and calculated cell volumes is 
less than 5 meV/H2. This error only has a minor effect 
on the relative total energies. 

The optimized B-H and Be-H bond lengths are given 
in Fig. [TJ As references, the B-H bond length in a 
(BH4)^ anion is 1.21 A, whereas a B-H bond length in 
a typical three center B-H-B bond is 1.34 A.'^'^ Compar- 
ison with these numbers indicates that the B~H bonding 
in Be(BH4)2 is closer to that in the (BH4)~, although 
there is some distortion due to the presence of the Be 
atom, in particular on the B^-Hc bond. This could in- 
dicate some competition between B and Be for bonding 
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FIG. 3: Total energy per unit cell of Be(BH^)2 as a function 
of the cell volume. 



TABLE I: Optimized atomic positions of Be(BH^)2. The 
space group is Micd (110) and all atoms are on Wyckoff po- 
sitions 16b. The optimized lattice parameters are a,b,c = 
14.28, 14.28,9.54 A. 



atom 


X 


y 


z 


Be 


0.2050 


0.0992 


0.0016 


Bd 


0.1695 


0.9702 


0.0068 


Bb 


0.1503 


0.1978 


0.1237 


Hd 


0.0997 


0.9439 


0.0653 


Hd 


0.2183 


0.9130 


0.9499 


H, 


0.2189 


0.0123 


0.0963 


H, 


0.1453 


0.0294 


0.9157 


Hi, 


0.1083 


0.1649 


0.0231 


Hi, 


0.2281 


0.1647 


0.1450 


Hi, 


0.1611 


0.2813 


0.1003 


Hi, 


0.1027 


0.1793 


0.2269 



to hydrogen. For comparison, the B-H bond lengths in 
alkah boranates are all very close to 1.21 A. The Be-H 
bond lengths are still quite large, however, which indi- 
cates a significant ionic contribution to the bonding. 

We have also optimized the structure of BeH2, see Ta- 
ble ini It agrees well with the experimental structure^^ 
and with that obtained in a previous calculation," the 
largest difference being that our calculated bulk modu- 
lus (21.4 GPa) is 10 % smaller than that calculated in 
Ref. [13 (23.8 GPa). For elemental boron we use the /3- 



TABLE II: Optimized crystal structure of BeH2. The space 
group is Ibam (72) and the optimized lattice parameters are 
a,b,c = 8.967,4.141,7.643 A. The experimental lattice pa- 
rameters are a,h,c = 9.082, 4.160, 7.707 k.^^ 



atom 


Wycoff 


X 


y 


z 


Be 


4a 








0.25 


Be 


8j 


0.1677 


0.1200 





H 


16k 


0.0882 


0.2241 


0.1520 


H 


8j 


0.3102 


0.2771 







FIG. 4: (color online) Charge density difference plots of 
Be(BH4)2 with respect to the isolated atoms. The top pic- 
ture gives a cut through a plane containing Bb, Hb and Be 
atoms of the polymer backbone, the bottom left picture a cut 
through the Bd, He side chain plane, and the bottom right 
picture a cut through the Bd, Hd side chain plane. 

rhombohedral structure as given in Ref. [3^ . For elemen- 
tal beryllium (space group P3ml (164)) we find lattice 
parameters a — 2.260 A and c = 3.567 A, which compare 
well to the experimental values of 2.29 A and 3.60 A, 
respectively!^ 



IV. ELECTRONIC STRUCTURE 

As discussed above the crystal structure of Be(BH4)2 
indicates a weak bonding between polymer chains, and 
a stronger bonding within a polymer chain. The charge 
displacement upon bond formation can be visualized by 
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plotting the charge density difference, i.e. the charge 
density of Be(BH4)2 minus that of the individual iso- 
lated atoms. Cuts through the charge density difference 
in various planes along a polymer backbone are shown 
in Fig. [31 They clearly indicate the formation of B-H 
covalent bonds, which are polarized somewhat towards 
the H atoms. The character of the Be-H bonds is much 
less clear from these plots. In any case these bonds are 
strongly polarized in the direction of the H atoms. 

The electronic projected density of states (PDOS) 
of Be(BH4)2, projected on s,p components of the in- 
dividual atoms, is shown in Fig. O Tetrahedrally 
bonded BH4 generates a characteristic pattern in the 
valence band part of the PDOS, which is qualitatively 
similar to that observed for AlH^ tetrahedra in the 
alanatesi^iiaiiiii^ii^iii The tetrahedral geometry of BH^ 
results in a splitting into two valence peaks, the lower 
one of s (Ai) symmetry and the upper one of p (T2) 
symmetry, with a relative weight ratio of 1:3. Projected 
on atomic states, the s-peak then has contributions from 
H s and B s orbitals, and the p-peak has contributions 
from H s and B p orbitals. The p-peak can be split due 
to symmetry breaking caused by the crystal field. This is 
clearly observed in the lowest two panels of Fig. [5l show- 
ing the PDOS on the Bd and Hd atoms with the s peak 
at ~ — 7 eV, and a p doublet around ~ — 1 eV. The split- 
ting between s- and p-peaks is large (~ 6 eV), and the 
crystal field splitting is much smaller 1 eV). 



The interaction between the BH4 units in the crys- 
tal lattice leads to a broadening of the peaks due to 
band formation. The interaction is strongest along the 
B;,(Hf,)4-Be-Bfc(Hh)4 polymer backbone, see Figs. [1] and 
This leads to an s-type band in the range ^ —9 to 
~ 7.5 eV, involving contributions from H;,, B^, and Be 
s orbitals, whose DOS has the characteristic shape of a 
one-dimensional structure, see the upper three panels of 
Fig. [5] In the range ~ — 5 to ~ — 2 eV we find and a 
set of p-type bands. The band widths are smaller than 
the sp splitting, but they are not negligible, reflecting the 
covalent bonding along the polymer backbone. 

The involvement of the Be atoms can be clarified by 
calculating the DOS for a Be(BH4)2 structure in which 
the Be atoms are replaced by a homogeneous background 
with charge 2+. The result is shown in Fig [51 The s 
and p valence bands discussed above disappear and are 
replaced by much narrower peaks that refiect electron 
localization on BII4 ions in this artificial structure. In 
other words, the Be atoms in Be(BH4)2 are involved in 
the covalent bonding. This is in contrast to alkali or 
alkaline earth boranates and alanates, where the DOS 
changes little if the cations are replaced by a background 
charge. The bonding in the latter compounds can be 
described as an ionic bonding between BH^ or AlH^ 
anions and M+ (akali) or M'^+ (alkaline earth) cations4i 
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FIG. 5: (color online) The electronic projected densities of 
states (PDOS) of Be(BH^)2. The Fermi level, Ef, at the top 
of the valence band is the zero of energy. The upper (red) 
and lower curves (black) give projections on p and s atomic 
states, respectively. Atomic radii of 0.7, 0.5 and 1.1 A are 
used for Be, B and H. 
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FIG. 6: (color online) The total electronic densities of states 
(DOS) of Be(BH4)2. The Fermi level, Ep, at the top of the 
valence band is the zero of energy. In the lower curve the 
beryllium atoms are replaced by an homogeneous background 
with charge 2+. 
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V. REACTION ENTHALPIES 

For light elements such as hydrogen, beryllium and 
boron the quantum character of their atomic vibrations 
is important. This leads to vibrational energies at zero 
temperature that are not negligable. For each compound 
involved in the reaction we calculate its zero point vibra- 
tional energy (ZPVE) from the frequencies of the vibra- 
tional modes in the optimized structure. For hydrogen 
molecules also the zero point rotational energy (ZPRE) 
is not entirely negligable. Reaction enthalpies AiJ at 
T = OK are then calculated from 

P r 

(1) 

where £^p°* denotes the total energy of the reaction 
products p or reactants r, E'^J^^ are the correspond- 
ing ZPVEs, and E'ff is the ZPRE of the hydrogen 
molecules involved in the reaction. 

For the hydrogen molecules we calculate a vibrational 
frequency of 4356 cm~^, in good agreement with the ex- 
perimental value of 4401 cm^^.i^ The ZPVE, 0.266 eV, 
is then calculated from the energy levels of a Morse po- 
tential. 





1 





















where u is the vibration eigenfrequency and = 
4.57 eV is the dissociation energy. Assuming that ortho- 
and para-hydrogen are produced in a proportion of three 
to one, the average ZPRE of a hydrogen molecule is 
0.011 eV, using the energy levels given in Ref. |45|. 

The calculated total energies and ZPEs of all com- 
pounds involved in the reactions are listed in Table [1111 
We consider two possible reaction paths for the forma- 
tion of Be(BH4)2. In the first path Be (6114)2 directly 
formed from the elements. 

Bc + 2B + 4H2(.9) ^Be(BH4)2. (3) 

The second path involves the formation of an intermedi- 
ate compound BeH2. 

Be-|-H2(5) ^BeH2. (4) 

BeH2 + 2B + 3H2(,g) ^Bc(BH4)2- (5) 

The enthalpies of these reaction are calculated using 
Eq. Il]) and the values given in Table IIIII 

Equation ([3]) gives a reaction enthalpy of —0.39 eV/H2 
if ZPEs are neglected. If ZPEs are included the reaction 
enthalpy becomes —0.12 eV/H2, which indicated the im- 
portance of ZPE corrections for these lightweight com- 
pounds. In principle these values are in a range that 
is useful for hydrogen storage. Using the ionic model 
of Ref.^* gives a reaction enthalpy of —0.02 eV/Il2 (ne- 
glecting ZPEs). In the previous section we have already 



TABLE III: Total energies (with respect to non spin polarized 
model atoms), zero point vibrational energies (ZPVE) and 
zero point rotational energy (ZPRE) in eV/formuIa unit in 
the relaxed structures. 





ETOt 


-gZPVE 


gZPRE 


Be(BHj2 


-43.353 


2.450 




BeHa 


-10.797 


0.542 




Ha 


-6.803 


0.266 


0.011 


B 


-6.687 


0.126 




Be 


-3.729 


0.091 





concluded that the bonding in Be(BH4)2 is not purely 
ionic, i.e. Be^"'"(BIl4)^. The Be atoms are bonded (par- 
tially) covalently to BH4, which increases the bonding as 
compared to the pure ionic picture, resulting in a higher 
dehydrogenation enthalpy. 

VIost alanates and other boranates form a simple al- 
kali / alkaline earth hydride when hydrogen is released 
in a first step. The dehydrogenation of this simple hy- 
dride then occurs as a separate second step. Usually the 
enthalpies are such that only the first step is considered 
useful for hydrogen storage. For Be(BH4)2 these two 
steps correspond to the reverse reactions of Eqs. ([5]) and 
(ID). The calculated reaction enthalpies of Eqs. ^ and 
© are -0.27 and -0.43 eV/Ha without ZPEs, and -0.09 
and —0.13 eV/H2 with ZPE corrections, respectively. 

Note that per H2 Be(BH4)2 is slightly more stable than 
BeH2. This would indicate that a one-step reaction di- 
rectly from the elements, Eq. ([3]), is more favorable than 
the two-step reaction via the simple hydride, Eqs. ^ 
and ([5]). The enthalpy difference however is very small. 
In addition kinetic barriers may influence the relative im- 
portance of the two reaction paths. 

We will now focus on finite temperature properties. 
For the solids we calculate the Gibbs free energy G{T) in 
the harmonic approximation 

G{T)^E^°* +H'''^{T)-TS^'^{T) (6) 

with 

H^'^{T) = ^ dujg{uj) ji^w -I- hujn{uj)j (7) 

and 

poo 

S'''^{T) = kB dwg{oj) {(3nwn{oj) - In [l - e"^'*"] }(8) 

where g{uj) is the phonon density of states, n[ijj) — 
[exp(/3fiii;) — 1] ^ is the Bose-Einstein occupation num- 
ber and (3 — l/ksT. The first term in the integral of 
Eq. ([7]) gives the ZPVE and the second term gives the 
finite temperature contribution. Note that we neglect 
the PV term (i.e. the distinction between energy and 
enthalpy), which is a good approximation for solids. 

For the Gibbs free energy, the enthalpy and the en- 
tropy of the hydrogen gas, we use the values given in 
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FIG. 7: (color online) Gibbs free energy of Be(BH^)2 and the 
constituting elements. 



Ref. H^. Fig. [7] gives the free energies of Be(BH4)2 and 
the products of the dehydrogenation reaction, i.e. the 
left and right hand sides of Eq. ([3]), at the standard pres- 
sure of 1 bar. At 162 K the free energy of the products 
drops below that of the hydride phase. This intersection 
of the two free energy curves defines the decomposition 
temperature. 

The predicted desorption temperature and equilibrium 
pressures rely only on the thermodynamics of the re- 
action. From the fact that experimentally Be(BH4)2 
seems to be stable at room temperature one may con- 
clude that kinetic barriers play an important role in sta- 
bilizing Be(BH4)2. In this respect Be(BH^)2 is similar to 
other boranates and alanates, where decomposition tem- 
peratures are much higher than what is expected on the 
basis of thermodynamics and catalysts have to be applied 
in order to overcome kinetic barriers. 



covalent as well as ionic contributions. The crystal struc- 
ture and the electronic density of states give evidence for 
[-BH4-Be-]„ (hehcal) polymers. 

The enthalpies of possible formation reactions are cal- 
culated including zero point energy corrections. The lat- 
ter are obtained by the calculating the phonon frequen- 
cies of all compounds involved in the reactions. Since 
not only hydrogen, but also boron and beryllium are rel- 
atively light elements, the zero points energies are rela- 
tively large for these compounds. The enthalpy of for- 
mation of Bc(BH4)2 from the elements is —0.39 eV/H2 
without and —0.12 eV/H2 with zero point energy contri- 
butions. 

Be(BH4)2 differs from other boranates and alanates in 
that its dehydrogenation to the elements is thcrmody- 
namically slightly more favorable than dehydrogenation 
via the simple hydride BeH2. In alkali or alkaline earth 
boranates and alanates dehydrogenation always occurs 
via the alkali or alkaline earth simple hydride. The dif- 
ferent behavior of Bc(BH^)2 is mainly caused by the high 
stability of bulk beryllium metal. 

Be(BH4)2 follows the general trends in the formation 
energies that have been observed in alkali and alkaline 
earth alanates and boranates. Boranates are more sta- 
ble than the corresponding alanates, lighter cations give 
compounds that are more unstable, and alkaline earth 
compounds are more unstable than alkali compounds.— 
Indeed Be(BH4)2 is less stable than LiBH4 or Mg(BH4)2. 
We have not found mentioning of beryllium alanate in 
the literature, which might indicate that this compound 
would be too unstable. 

Using the calculated phonon spectrum we have calcu- 
lated free energies within the harmonic approximation to 
assess thermodynamic properties at finite temperature. 
We obtain a decomposition temperature of 162 K at a 
1 bar pressure. 



VI. CONCLUSIONS 
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